Molecular spring: from spider silk to silkworm silk 
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In this letter, we adopt a new approach combining theoretical modeling with silk stretching 
measurements to explore the mystery of the structures between silkworm and spider silks, leading to 
the differences in mechanical response against stretching. Hereby the typical stress-strain profiles are 
reproduced by implementing the newly discovered and verified "/3-sheet splitting" mechanism, which 
primarily varies the secondary structure of protein macromolecules; Our modeling and simulation 
results show good accordance with the experimental measurements. Hence, it can be concluded 
that the post-yielding mechanical behaviors of both kinds of silks are resulted from the splitting of 
crystallines while the high extensibility of spider dragline is attributed to the tiny /3-sheets solely 
existed in spider silk fibrils. This research reveals for the first time the structural factors leading to 
the significant difference between spider and silkworm silks in mechanical response to the stretching 
force. Additionally, the combination of theoretical modeling with experiments opens up a completely 
new approach in resolving conformation of various biomacromolecules. 

PACS numbers: 87.85. J-, 81.70.Bt, 87.15.ap 



Natural silk spinning was pioneer of economy flourish- 
ing worldwide since the first industrial revolution, and 
was continuing acting as a pillar industry during the last 
several centuries until massive manufacturing of man- 
made fibres derived from feedstocks of petrochemical was 
popularized. Spider silk, superior to other biomateri- 
als, has extraordinary strength comparable to steel, and 
the highest toughness among all natural silk fibres up to 
date [J 0, d, 0]. It is suggested that one strand of pencil 
thick spider silk can stop a Boeing 747 in flight. The ex- 
ceptional properties serve new application in industry to 
fulfill various functions, such as bullet-proof vests, rein- 
forced composites and aircrafts panels substitute. On the 
other hand, the high production of silkworm silk makes 
its irreplaceable stay in textile and other markets, despite 
its inferior properties to spider silk[B|, [||. Complementar- 
ily, the sustainable application has been penetrated in 
daily life, ranging from costume manufacture to clinical 
treatment. Both two kinds of silks exhibit greater envi- 
ronmental friendliness and bio-compatibility than man- 
made petrochemical materials, which implies the conve- 
nience of its synthesis, fabrication and recycle Al- 
though spider and silkworm silk share a high degree of 
similarity in their chemical composition and microscopic 
structures, the two kinds of silks behave differently in 
their mechanical responses 0, [1(1 EH S3]- To fulfill their 
potential functionalities requires a full understanding of 
the underlying mechanism; however, to date very few 
works have been engaged in this exploration. 

Both spider silk and silkworm silk are primarily com- 
prised of Alanine- and Glycine-rich polypeptides in the 
fibrils; Besides, natural spun silkworm silk contains two 
strands coated with sericin, and the spider dragline 
silk contains major ampullate(MA) and minor ampul- 



late(MI) silk. X-ray diffraction and AFM probing[l3| 
provide direct evidence on the molecular level that subtle 
deviation in amino acids assembly sequence may lead to 
significant difference in mechanical properties [lj, 15, lij ]. 
However, it is still not sufficient to fully explain the me- 
chanical difference. As can be seen, the stress-strain pro- 
file of silkworm silk is segmented in two regions, sep- 
arated by yielding point. In low stress regime, the 
silk fibre behaves linearly elastic and skips into the 
post-yield re gion with glass state which is nonlinear and 
irrecoverable [18j. Different from silkworm silk, spider 
draglines show the work-hardening phenomenon in the 
post-yield region, describing its dramatic increase of the 
elasticity of spider silk when subjected to certain extent 
of stretch, after which the fibres turn softer again. 

Accordingly, the mechanical differences of spider 
draglines and silkworm silk are also related to the high- 
order structures of protein macromecules (R [rj [TBI 1^ . 20j |. 
The key problem is to determine structural factors, which 
affect mechanical response to the external stress and 
cause the differences of the two kinds of silks. The 
secondary structure of silks can be divided into crys- 
talline, primarily the /3-crystallites, and non-crystalline 
(amorphous) domains (Fig. 1(a), (b)), which contains 
random coils, a-helices and /3-sheets |2ll. [22j . It is be- 
lieved that the secondary structures are similar in both 
two kinds of silk fibres except that ordered /3-sheets were 
solely observed only in the amorphous region of spider 
draglines 13. |20|. [23|. certainly, the crystallinity, denot- 
ing fraction of crystalline region, and the sizes of (3- 
crystallites are different. We notice that one /3-sheet 
in the amorphous region of spider draglines are stacked 
by part of a single protein macromolcculc, while the (3- 
crystallites in the crystalline regions are formed by sev- 
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FIG. 1: The left panels are the structures of fibrils aligning 
along silkworm silk and spider silk fibres in (a) and (b) re- 
spectively. The upper-right sketches are the nano-structures 
of each fibril; the larger blocks denote the inter-protein f3- 
crystallites, and the smaller beads denote the intra-protein j3- 
sheets; the flexual lines in-between denote helical or random 
structures in the amorphous matrix. The lower-right parts is 
the model construction based on their nano-structures, and 
the element correspondence is directed by the thin solid ar- 
rows. 



eral adjacent of them. In order to distinguish the two 
kinds of /3-sheet structures, the one in crystalline region 
is referred to the "inter-protein (3 - crystallite" , and the 
one in amorphous region is referred to the 11 intra-protein 
/3 - sheet" , as shown in Fig. 1(a) and (b). 

In this letter, aiming at revealing the connection be- 
tween the mechanical properties of spider and silkworm 
silk and their nano structures, we are inclined to un- 
derstand the roles of the elements based on our models. 
Therefore, we are observing the collective behaviors of 
two kinds of silk, in terms of the structural difference, 
typically for the stress-strain profile. 

Let us first start our modeling of silkworm silk. As 
shown in Fig. 1(a), /3-crystallites are extracted as massive 
bodies, and the in-between springs represent the amor- 
phous matrix, including random coils, a-helices, etc. A 
quasi-periodic structure(Figure la) is adopted with re- 
peated segments of rigid bodies and springs, and they 
are aligned along fibre axis subjected to the one dimen- 
sional confinement we simplified to. As indicated by Ra- 
man spectrum experiments, the /3-crystallites and amor- 



phous regions are more likely to connect in series ll| , 
which is well described by the quasi-periodic structure. 
In general, the transition from rubber to glass state at 
the yielding point requires well definition in both the two 
regions. Hence, Hooke's springs are imposed to describe 
the linear elastic behavior contributed mostly by amor- 
phous region before reaching the yielding point, and be- 
come elasto-plastic with nonlinearity and irreversibility 
in the post-yield region. In addition, /3-crystallites, pre- 
viously treated as rigid bodies, are now actually elas- 
tic with relatively 4 times high elasticity of the whole 
fibre 0. Taking into account the facts of the transverse 
size of /3-crystallites shrinking and the irreversible pro- 
cess of elasticity measurement for the pre-stretched silk 
fibres, "/3-sheet splitting" turns to be possible mechanism 
in characterizing the yielding behavior 1^, 11 | , and will 
be discussed in the following context. 

Based on aforementioned facts, the model can be ex- 
plicitly described as a serially connected system of N 
segments, each of which is comprised of massive body 
with mass m and light spring with original elasticity ko. 
The segment at one end is stretched by an applied force 
F and the other boundary one is fixed immobile. During 
the stretch process, for an arbitrary segment i, the exten- 
sion linearly increases with the responsive force F(i) in 
the linear region until F(i) reaches a criteria of threshold 
force F t h{i) , after which the /3-crystallites start to split. 
The collective behavior of splitting at {F t h(i)} charac- 
terizes the critical behavior around the yielding point. 
Because the splitting of /3-crystallites requires external 
energy input equal to or greater than the cohesive energy 
of hydrogen bonds, equivalently a large applied force F 
is likely to break more hydrogen bonds, thus to retrieve 
longer amorphous protein molecules from the compactly 
stacked /3-crystallites. Hereby, a simple linear relation 
characterizing such observation between applied force F 
and the released molecule length AL(i) can be presumed 
as: 



AL(i) 



F < F th (i) 
F th (i))/E F>F th (i) 



(1) 



Therefore, the relation between applied force F and 
observable extension Ax (— £i=i & x i) carl be deduced 
in Eq.(2) as: 
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1 



-AL) 
(Ax- 



il) 



= 1 kt 



N 

E 

i=l 



AL(i)) 



= k E L Ax - Nk L ■ f[ (F - f) Pth (f)df 
NE L Q + Nf F (F-f)p th (f)df 

where pth(-) is the distribution of the threshold force 
{Fth(i)}, Lq is the original length of each segment, and 
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FIG. 2: The relation between stress a and strain e of silkworm 
silk fibres: the black solid line is the experimental curve, and 
the red dashed line is the simulation results. A computer- 
controlled motorized spindle was used to draw fibers from 
silkworm with 12mm/sec, 24mm/sec and 36mm/sec from 
the bottom up. Two successive 45mm, 95°C heating proce- 
dures in 0.5% Na,2COz and 1% soap solution was adopted in 
the degumming process. Measurements were performed using 
an Instron MicroTester (Model 5848; force resolution, 0.5% 
of indicated load; position resolution 0.02/im; strain rate is 
50%/min), at 20±2°C and the humidity was kept at 60±5%. 
the scaling constant in simulation S c — 20, the corresponding 
parameters ko, Eq, and (F t h) increase with the reeling speed. 



h(x)(J^ oo 5(x)dx) is the step function. Besides, the cal- 
culation of elasticity k c s follows the law of sequentially 
linked springs with effective elasticity of each segment 
ki = ko/(l + AL(i)/ Lq). Additionally, for large system 
size N, continuity approximation is applied in deriving 
Eq.2. Mapping three dimensional reality to our one di- 
mensional model, the scaling is governed by such regu- 
lations: F ~ er(set the area of cross section AS = 1), 
and Ax ~ e • NLq, where a and e denote the tensile 
stress and tensile strain respectively. As one of the key 
features of biomaterial, ununiformity caused by random 
generated defects leads to serious instability of their me- 
chanical properties. Without loss of generality, gaussian 
random function is chosen for the distribution of thresh- 
old forces. 

To mimic the dynamic process of silk stretching, the 
extensive numerical simulations based on Molecular Dy- 
namics (MD) were carried out. In comparison, the 
stress-strain profiles of silkworm silk(Fig.2) and spider 
draglines(Fig.4) under different reeling speeds were mea- 
sured. As shown in Fig. 2, MD numerical results of stress- 
strain profiles coincide very well with experimental data 
at the proper values E /ko and (F t h)- It is found that the 
ratio Eo/ko determines the relative steepness of the two- 
segmented stress-strain profiles, with two limits that the 
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FIG. 3: (a):Stress-strain profile of silkworm silk for different 
spread AS of the gaussian form pth(-) and same expectation 
value (Fth)- (b):Stress-strain profile of silkworm silk for dif- 
ferent forms of pth(-) with same {Fth) and various spread. 



post-yield behavior is strain-independent for Eo/ko ~ 0, 
and remains linear elastic for Eo/ko ~ oo. The stress- 
strain profile is relatively robust to the form of distri- 
bution function of threshold forces, as general gradually 
varying and centralized functions give out similar profiles. 
As shown in Fig. 3, various distribution of threshold force 
Pth gives no essential difference, but only little variation. 
Additionally, different reeling speed corresponds to their 
specific (Fth) an d Eo/ko, reflecting the initial formation 
difference of the nano structures of macromolecules, pre- 
sumably related to the /^-crystallites position and orienta- 
tion. Therefore, the breaking of /3-crystallites in silkworm 
silk fibrils weakens the linkage between protein molecules, 
which leads to its softening in the post-yield regime. 

To model the structure of spider dragline, we take into 
account the performance of the intra-protein /3-sheets in 
the amorphous region (Fig. 1(b)), which is uniquely ob- 
served spider draglines, in addition to adopting the idea 
of the above model. Analogous to the silkworm silk mod- 
eling, the linear assumption of protein molecule length 
split from intra-protein /3-sheets {AL'(i)} versus stretch- 
ing force F is preserved, with proportional constant E\ 
and threshold forces of splitting {F{ h (i)}. Due to the 
morphological inperfectness and stacking incompactness 
of the tiny intra-protein /3-sheets, they are easier to be 
split than the inter-protein /3-crystallites, thus to release 
longer protein molecules given the same external energy 
consumption. Besides, the limited size of intra-protein j3- 
sheets lead to their complete destruction during stretch- 
ing at forces {Ft r (i)}. Hereby, if the stretching force 
F larger than this terminating force, the intra-protein 
/3-sheets can no longer release any length of protein 
molecules; instead, they are fully destroyed, and behave 
similarly as coiling structures in the amorphous matrix. 
Concerning the above statements, qualitative relation 
can be derived as E\ < E , and (F{. h ) < (F tr ) < (Fth), 
which should be well obeyed in the following modeling. 
Force-extension relation of spider draglines can be anal- 
ogously derived as that of silkworm silk, by setting the 



4 



similar form of the lengthening term {AL(i)} for inter- 
protein /3-crystallites in Eq.l, and a relatively more com- 
plicated form for intra-protein /3-sheets shown in Eq.3. 

f F < F[ h (i) 

AL(i) = UF- Fl h (i))/Ex Fl h (i) <F< F tr {i) 
{{F tr {i)-Fl h {i))/E l F>F tr (i) 

(3) 

The explicit force-extension expression of spider silk 
turns out to be: 

_ kpLpAx - Nk L ■ Q(F) 

NL + NCl(F) ' 1 ' 

where f2(F) denotes: 

= ^(i- j\tr{fWj £ (F - f)p' th (f)df 

+ ^ [ F (F-f) Pth (fW- (5) 
Eo Jo 

Therein, p' th {-), ptr{-) and pth(-) respectively charac- 
terize the distribution of the threshold forces, the criti- 
cal forces terminating the intra-protein /3-sheets splitting, 
and the threshold forces for inter-protein /3-crystallites. 
Correspondingly, Eq.5 shows the contribution from intra- 
protein /3-sheets splitting, their complete destruction and 
the splitting of inter-protein /3-crystallites. The occur- 
rence of each mechanism is also reflected in the segmenta- 
tion of stress-strain profile for spider draglines. Similarly 
as silkworm silk, simulations have been carried out to re- 
produce the stretch process of spider draglines. As shown 
in Fig. 4, the parameter ratio E\/(pkp) and Eq/((1— p)kp) 
are dominant to characterize the relative steepness of 
stress-stain profile. Besides, the stress-strain profile is 
robust to the exact form of force distribution functions; 
instead, only their central values denoting general each 
mechanism in effect need to be determined. Moreover, it 
follows that linear elasticity kp and proportional constant 
Ep and E\ increases with reeling speed, indicating their 
dependence of initial conformation of structures. There- 
fore, the mechanism of each segment of the stress-strain 
profile of spider dragline has been clearly distinguished: 
the intra-protein /3-sheets split at the first yielding, and 
their complete destruction gives rise to the lengthening of 
the fibrils and the stiffening in the post-yield region, thus 
the work-hardening phenomenon is well reproduced. On 
the other hand, the inter-protein /3-crystallites split dur- 
ing or after the release of intra-protein /3-sheets in the 
similar way of silkworm silk, resulting in their soften- 
ing before breaking. The diverse high-order structures of 
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FIG. 4: The black solid lines denote the stress-strain profile 
for spider silk under different reeling speed [3]. Stress-strain 
curves of spider dragline silk (2.5, 10,25, lOOmm/sec motor- 
reeled at 22° C, from the bottom up) were performed using 
an Instron MicroTester (Model 5848; force resolution, 0.5% 
of indicated load; position resolution 0.02/im; strain rate is 
50%/min), at 22°C and the humidity was kept at 55 — 60%. 
The red dashed lines are the computational results, where 
scaling constant S c = 20, the preserved fraction of intra- 
molecule /3-sheet p — 0.5, and parameters ko, Eo and Ei 
are varied for different reeling speeds. 

protein molecules for silk fibres of different species result 
into their distinct mechanical complexity. 

In conclusion, we have established the correlation be- 
tween the nano structures and mechanical properties of 
silk fibres in our essay, on the basis of modeling in combi- 
nation with the silk stretching experiments. We obtained 
for the first time the breaking mechanism of the two 
types of silks. The splitting of inter-protein /3-crystallites 
gives rise to the weakening the linkage among molecules 
while limited capacity of intra-protein /3-sheets results 
in the lengthening and the extra stiffening in the post- 
yield regime of spider draglines. Moreover, the explo- 
ration from molecular mechanical properties turns out 
be a brand new probing method to investigate the sub- 
tle structural differences among biomacromolecules. This 
provides a new passage in understanding the origin of 
molecular mechanical behaviours, and will facilitate the 
identification of robust technologies in fabricating silks of 
ultrafunctionality. 
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